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A model for three-dimensional Rayleigh-Benard convection in low-Prandtl-number fluids near on¬ 
set with rigid horizontal boundaries in the presence of a uniform vertical magnetic field is constructed 
and analyzed in detail. The kinetic energy K, the convective entropy $ and the convective heat flux 
{Nu — 1) show scaling behaviour with e = r — 1 near onset of convection, where r is the reduced 
Rayleigh number. The model is also used to investigate various magneto-convective structures close 
to the onset. Straight rolls, which appear at the primary instability, become unstable with increase 
in r and bifurcate to three-dimensional structures. The straight rolls become periodically varying 
wavy rolls or quasiperiodically varying structures in time with increase in r depending on the values 
of Prandtl number Pr. They become irregular in time, with increase in r. These standing wave 
solutions bifurcate first to periodic and then quasiperiodic traveling wave solutions, as r is raised 
further. The variations of the critical Rayleigh number Raoa and the frequency cuoa at the onset of 
the secondary instability with Pr are also studied for different values of Chandrasekhar’s number 

Q. 

PACS numbers: 47.35.Tv, 47.20.Bp, 47.20.Ky 


I. INTRODUCTION 

The role of magnetic field on thermal convection in 
low-Prandtl-number fluids [im3 has received consider¬ 
able attention because of its importance in geophysical 
and astrophysical problems. The vertical magnetic field 
delays the onset of primary instability [11,9 , while the 
horizontal magnetic field dela ys t he onset of oscillatory 
(secondary) instability H, 9 Till . The latter also 
makes the straight rolls align along its direction. The 
external magnetic held also reduces the convective heat 
flux 9 [l9 Ill, [13, [13|- The scaling behaviour of ki¬ 
netic energy, convective entropy and heat flux near on¬ 
set were studied recently with stress-free top and bottom 
plates [131. They are not investigated for the magne¬ 
toconvection with rigid horizontal plates. The secondary 
instability also leads to interesting three dimensional dis¬ 
sipative structures. It is assumed that the Nusselt num¬ 
ber Nu and the mean convective entropy $ = ^ / 9‘^dV 
also scale like the mean kinetic energy K = ^ f v^dV. Re¬ 
cent numerical simulations [13|, however, show that Nu 
scales with e = [Ra/RadQ) — 1] as e“ with a < 1 near 
the onset, where Ruc is the critical value of the Rayleigh 
number at the primary convection. 

In this paper, we present a model of Rayleigh-Benard 
convection in metallic fluids in the presence of uniform 
vertical magnetic field with rigid, thermally conducting 
and electrically insulating horizontal plates. We use the 
model to investigate various properties of magnetocon¬ 
vection near onset. The kinetic energy K is found to 
be proportional to e near onset. The convective entropy 
<I> and the convective heat flux (Nu — 1), scale with e 
as ® and respectively. The oscillatory convec¬ 

tion appears at the secondary instability as r is raised 
slowly. The distance from the criticality (roe — 1) for os¬ 


cillatory instability scales with the Prandtl number Pr 
as Pr^-^. We have observed the possibility of traveling 
as well as standing wave solutions above the onset of sec¬ 
ondary instability. They were found to be periodic as 
well as quasiperiodic in time. A new traveling convec¬ 
tive pattern consisting of alternating oblique wavy rolls 
is also observed. The time averaged root mean square of 
the velocity and temperature fields are also investigated. 


II. HYDROMAGNETIC SYSTEM 

We consider a thin horizontal layer of electrically con¬ 
ducting Boussinesq liquid of mean density po, thickness d, 
kinematic viscosity v, thermal diffusivity k, magnetic dif- 
fusivity A, and thermal expansion coefficient a confined 
between two parallel plates, which is uniformly heated 
from below and uniformly cooled from the top. An ad¬ 
verse temperature gradient /3 is maintained across the 
fluid layer and a uniform magnetic field Bq is applied 
anti-parallel to the acceleration due to gravity . We have 
chosen a Cartesian coordinate system with origin in the 
middle of the fluid layer, the xy- plane coincident with 
the horizontal plane. The unit vector es is directed along 
the vertically upward direction, which is considered to be 
the positive direction of the z- axis. The fluid is initially 
at rest and allows conduction of heat flux along the verti¬ 
cal direction. As soon as /? is raised above a critical value 
Pc keeping Bq fixed, magnetoconvection sets in. We ig¬ 
nore the magnetic Prandtl number Pm = irjX which is 
generally of the order of 10“^ or smaller for terrestrial 
fluids. The magnetohydrodynamics is then governed by 
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the following dimensionless equations: 


dtv + {v ■ S/)v = — Vp + + Qdzb + RaOes, (1) 

= -dzv, (2) 

Pr[dte + {vW)0]=V^0 + V3, (3) 

V • u = V • 6 = 0, (4) 

where v {x,y,z,t) = ivi,V 2 ,V 3 ) is the velocity field, b 
{x,y,z,t) = (&i, 62 ,^ 3 ) is the induced magnetic field due 
to convection, 0 (x, y, z, t) is the deviation in the temper¬ 
ature field from the steady conduction profile. Lengths, 
time, temperature field and induced magnetic field are 
measured in units of fluid depth d, viscous diffusion time 
vlid/n and Bqu/X, respectively. The magneto¬ 
convection is controlled by the three dimensionless pa¬ 
rameters: (i) the Rayleigh number Ra = ajdgd'^ /vk, a 
measure of the buoyancy force, (ii) the Prandtl num¬ 
ber Pr = vjK, and (iii) Chandrasekhar’s number Q = 
B^d^/A'KpovX, which is a measure of the imposed mag¬ 
netic field. The rigid, thermally conducting and elec¬ 
trically insulating horizontal plates leads to following 
boundary conditions [ 2 ^: 

Bb 

vi=V 2 = V3=6 = 0, bi =b2 = = 0 (5) 

oz 

at z = ± 1 / 2 . 


V3{x,y,z,t) = [ITioi cos (fca:) ± ITioi sin (fcx) 

± Will {t) cos (kx) cos (ky) 

+ Will {t) sin (fcx) cos (ky) 

+ Wiii{t)cos{kx)sm{ky) 

+ lTni(0 sin (fcx) sin (fcy)]C'i(Aiz), ( 8 ) 


W 3 ix,y,z,t) 


= [Zoii{t)cos{ky) 

± Zqii {t) sin {ky)] cos(7rz) 

± [^112 {t) cos {kx) cos (ky) 

+ Z 112 {t) sin {kx) cos {ky) 

± Z 112 {t) cos {kx) sin {ky) 

+ Zii 2 {t) sin {kx) sin {ky)] sin (27rz), (9) 


0 {x,y,z,t) 


= [Tioi cos {kx) + TiQi sin {kx) 

+ Till {t) cos {kx) cos {ky) 

+ Tiii{t) sin {kx) cos {ky) 

+ Tiii{t) cos {kx) sin {ky) 

+ Tni(t) sin (fccc) sin (fcy)] cos (ttz) 

± roo2(<)sin(27rz), (10) 


III. THE MODEL 


where 


We now construct a low-dimensional model to investi¬ 
gate the essential features of Rayleigh-Benard magneto¬ 
convection with realistic boundary conditions near onset. 
First we eliminate the pressure term from Eq.[T]by taking 
curl once on both sides of Eq. [TJ The projection of the 
resulting equation on vertical axis gives an equation for 
the vertical vorticity W 3 , which is given as: 

dtUJ3 + {v ■ V)w3-(w • V)u3 = V^a;3±(5e3-[Vx(5^6)] (6) 

Operating by curl twice on Eq. [TJ using the equation 
of continuity, and then projecting on vertical axis leads 
to 


dtV'^V 3 ± 63 • Vx [(w • V)r; — (u • V)a;] 

= V‘^V 3 ± RaV]j0 - Qe 3 -[V x Vx(5^6)] (7) 

We now follow the procedure used by Niederlander et 
al. for making a model with no-slip conditions in the 
context of ferro-fluids. All the fields are taken to be peri¬ 
odic in horizontal plane. Chandrasekhar’s functions are 
used in the vertical direction for the expansion of ver¬ 
tical velocity field so that the no-slip conditions can be 
applied on velocity field at z = ± 1 / 2 . As the magnetic 
field is slaved to the velocity field, it can be easily com¬ 
puted using Eq.|2|once we know the velocity modes. The 
expansion for the vertical velocity U 3 , the vertical vor¬ 
ticity u;3 = (V X v)3 and the convective temperature 0 , 
compatible with the boundary conditions (Eq. |5|) are: 


Cl (Aiz) 


coshAiz cosAiz 
coshAi /2 cosAi /2 


( 11 ) 


is the first order Chandrasekhar’s function d, [ 2 ^ with 
Al ~ 4.73. The horizontal velocities are then computed 
using the relations: 


= -dxzV3 - dyU }3 (12) 

= -dyzVs + dxUJ 3 . (13) 

Once all the three components of the velocity field are 
known, the horizontal vorticities can be computed easily. 

By projecting the hydromagnetic system of equations 
[Eqs. [T]|3| on these modes, we get a model for magneto¬ 
convection with no-slip, thermally conducting and elec¬ 
trically insulating boundary conditions. The model con¬ 
sists of nineteen coupled ordinary differential equations. 
The shear flow, in general, can not be fully expressed in 
terms of the vertical velocity and the vertical vorticity 
[ 1 ^ [30l| . For example, the dimensionless shear stresses 
o'la = ± dsvi and (T 23 = 82 X 3 ± 83 X 2 are nonzero, 

even if a part of the horizontal velocities vi and V 2 depend 
only on the vertical (the third) coordinate. However, the 
terms dependent purely on the vertical coordinate cannot 
be generated from the vertical vorticity W 3 . The model 
presented here considers the possibility of the shear gen¬ 
erated by the vertical vorticity but ignores the generation 
shear due to the perturbations independent of W 3 . 
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The model is integrated by the standard fourth or¬ 
der Runge-Kutta (RK4) method with a dimensionless 
time step of 10“^. We set k = kc{Q), where the criti¬ 
cal wave number kc{Q) is known from Chandrasekhar’s 
linear theory Q . We first determine the critical Rayleigh 
number Rqc from the model and compare them with the 
well known results of Chandrasekhar. Table |T] enlists the 
values of critical Rayleigh number RadQ) for different 
values of Q as obtained from the model with the val¬ 
ues known from the linear theory Q. For (5 = 0, the 
value of Rac is 1728, which is exactly equal to the value 
obtained by Niederlander. This value is within 1.2% of 
the value obtained by the linear theory Q. The error in 
determination of the critical Rayleigh number from the 
model is within 5% for Q = 50. We restrict ourselves 
upto Q = 100 for which the maximum error in Rac is 
less than 6.6%. 


Q 

kc{Q) 

Linear theory 
RadQ) 

Model 

RadQ) 

Error 

0 

3.13 

1707.8 

1728 

1.18% 

10 

3.25 

1945.9 

1926 

1.02% 

50 

3.68 

2802.1 

2664 

4.93% 

100 

4.00 

3757.4 

3510 

6.58% 


TABLE I: A comparison of the critical Rayleigh numbers ob¬ 
tained from Chandrasekhar’s linear theory and the model. 


IV. STABILITY OF STRAIGHT ROLLS 


For a thin layer of metallic fluid {Pr > Pm) confined 
between rigid, thermally conducting and electrically in¬ 
sulating horizontal boundaries, the magnetoconvection 
appear as stationary straight rolls, as in the case of free- 
slip boundaries. The vertical magnetic field delays the 
onset of stationary convection. We now find the fixed 
points and their stability by analyzing the model. Just 
above the onset of straight rolls there are only five non¬ 
zero modes: LFioi, kFioii ^loi, ^loi a-^d 0oo2- The roll 
fixed points are given by the following relations: 


= /i(c Q, Pr)ei^i, TTfoi = h{r, Q, Pr)e*,„ 

0*002 = f2ir,Q,PT), 011^+0*1^ = -^, (14) 


where 


/i(r,Q,Pr) 

f2{r,Q,Pr) 


_ airRac{Q)kl{Q) _ 

0,2 + a-ikl{Q) + aik‘^{Q) + a^Q 
-hi + b2h-hkl{Q) 
bifiPr 


(15) 


with ai = 1.376 x 10i°, az = 9.881 x ag = 4.857 x 
10“, 04 = 1.974 X 10^°, 05 = 1.895 x 10“, bi = 9.741 x 
101°, b2 = 1.377 X 101°, 63 = 9.869 x 10°, 64 = 5.032 x 
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FIG. 1: (a) Real and (b) Imaginary parts of the largest eigen¬ 
value as a function of r for different values of Q {Pr = 0.1). 
The onset of oscillatory instability is delayed and the fre¬ 
quency at the onset increases with increase in the Chan¬ 
drasekhar’s number Q. 


10i°, c = 15.485. The temperature mode roo 2 is the only 
nonlinear mode for the fixed point in the form of straight 
rolls. We now compare these modes with the expansion 
of the fields given by Clever and Busse [l^ with no-slip 
horizontal boundaries. These five modes are exactly the 
same as the modes obtained by truncating the expansion 
of the velocity field at the first term and that of the 
convective temperature at the second term. Small errors 
in the critical Rayleigh number Rqc for small values of Q 
suggest that the model describes the qualitative features 
of magnetoconvection near the onset well. 

Rolls could be in any direction in an extended layer of 
metallic fluid. As the reduced Rayleigh number r is raised 
slowly in steps, the perturbations in the form of vertical 
vorticity may be excited through nonlinear interaction 
with the vertical velocity. We investigate the stability of 
rolls in the presence of additional fourteen modes. We 
find the eigenvalues of a 19 x 19 matrix, obtained by 
linearizing about the roll fixed points. The eigenvalue 
A„i with the largest real part is found to form a complex 
conjugate pair. 

Figure [T] shows the variation of (a) the real and (b) 
imaginary parts of Am with r for different values of Q 
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FIG. 2: (a) Real and (b) Imaginary parts of the largest eigen¬ 
value as a function of r for different values of Pr {Q = 10). 
The onset of oscillatory convection is delayed and the fre¬ 
quency at the onset decreases with increase in the Prandtl 
number Pr. 


for Pr = 0.1. The real part of Xm becomes positive at 
greater value of r for larger Q values. This suggests that 
the oscillatory instability via forward Hopf bifurcation is 
also delayed by the vertical magnetic held. The frequency 
ujos at the onset of oscillatory instability shifts to higher 
values, as Q increases. 

The onset of oscillatory instability is greatly affected 
by the variation of the Prandtl number Pr. Figure [2] 
displays the variation of (a) the real and (b) imaginary 
parts of Xm with r for different values of Pr for Q = 
10. The onset of oscillatory instability shifts to a larger 
value with Pr and the frequency at the onset decreases 
with increasing values of Pr, which is consistent with the 
behaviour in the absence of any magnetic Held. Figure [3] 
(a) displays the variation of the threshold for oscillatory 
instability Cos = fosiQ, Pt) — 1 with Prandtl number Pr 
for different values of Q. The threshold rgg scales with 
Pr as Pr^-^ for Pr > 0.1 and the scaling appears to be 
independent of Q. The value of ros is larger for higher 
values of Q at a Hxed value of Pr. The frequency to os at 
the onset of oscillatory instability increases with increase 
in Q at a Hxed value of Pr. The frequency uios decreases 
with Pr for Pr > 0.1. However, the variation of the 


FIG. 3: Variation of (a) the critical value of the reduced 
Rayleigh number ros{Q, Pr) —1 and (b) frequency at the onset 
of oscillatory magnetoconvection as a function of as a function 
of the Prandtl number Pr for different values of Q. 


frequency with Pr is non-monotonic for much smaller 
values {Pr < 0.1). 


V. SCALING OF GLOBAL QUANTITIES NEAR 
ONSET 

We then use the model to investigate possible scal¬ 
ing behaviour near the onset of magnetoconvection. We 
start integration of the hydromagnetic system with ran¬ 
domly chosen initial conditions for a given value of Pr 
and r. The value of r is raised in small steps keeping 
Pr Hxed. The Hnal values of all the Helds of the last 
run are then used as initial conditions of a fresh run. As 
Pr is always greater than Pm, which is assumed to be 
vanishingly small, we always observe stationary straight 
(two-dimensional) rolls at the primary instability. This 
is consistent with the Chandrasekhar’s prediction. We 
Hnd various time dependent dissipative structures at sec¬ 
ondary and higher order instabilities, as r is raised fur¬ 
ther. The whole process is repeated for different values 
of Pr. We have also checked several integration with 
random initial conditions to Hnd out any possibility of 
hysteresis. We did not find any hysteresis in the model. 
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FIG. 4: Variations of (a) the mean kinetic energy per unit 
mass K and (b) the mean convective entropy with the pa¬ 
rameter e = [Ra/Rac{Q) — 1] for Pr = 0.025, as obtained 
from the model. 


We now present the results of scaling behaviour of 
three global quantities: the kinetic energy per unit mass 
K = ^ f v^dV, the convective entropy per unit mass 
<& = i J 9^dV, and the convective heat flux {Nu — 1) 
with e = [Ra/RadQ) — 1] for different values of Chan¬ 
drasekhar’s number Q. The parameter e is a measure 
of the distance from the criticality. Figure IH^a) shows 
the the kinetic energy it' as a function e for Pr = 0.025 
and for four different values of Q. The kinetic energy 
increases linearly with e for time independent magneto¬ 
convection. This means that the average speed of the 
fluid flow is proportional to y/e at the primary instabil¬ 
ity, which is a well known result. The value of K is 
higher for larger values of Q. The kinetic energy shows 
a sharp decrease, as e is raised in small steps. The sharp 
decrease in the kinetic energy is accompanied by a time 
dependent magnetoconvection. The convection is found 
to be quasiperiodic in time at the onset of secondary in¬ 
stability. The time-averaged value of the kinetic energy 
first decreases, attains a minimum, and then increases 
once again with increase in e. The sharp decrease in K 
occurs at higher values of e for larger values of Q. The dip 
in K is also shallower at higher values of Q. The larger 
vertical magnetic field delays the onset of secondary in¬ 
stability, which is time-dependent. As e is raised further, 
time averaged value of K is again found to be varying 


almost linearly with e for 1 < e < 10. 

Figure lU^b) shows the variation of the convective en¬ 
tropy $ with e for Pr = 0.025. The qualitative behaviour 
is similar to that of the kinetic energy but with some sig¬ 
nificant differences: $ scales with e as ® near the onset. 
The convective temperature field is then proportional to 
gO.45 jg different than the behaviour of the aver¬ 

age speed. This is due to the fact that a part of the 
available thermal energy is used to maintain a net ther¬ 
mal flux in the vertically upward direction. There is no 
net momentum flux in the vertical direction. The scal¬ 
ing behaviours of the average kinetic energy K is slightly 
different from that of the convective entropy <f>. The 
entropy also shows a sharp decrease at the onset of time- 
dependent (secondary) instability. For higher values of 
e (1 < e < 10), all the curves for $(e) have a common 
slope approximately equal to 0.4, which is much smaller 
than its value (~ 0.9) just above the onset of magneto¬ 
convection. 



e 


FIG. 5: Plot of convective heat flux {Nu — 1) versus e = 
[Ra — Rac{Q)\lRadjQ) for Pr — 0.025 as obtained from the 
model (M) and compared with the results of Glever and Busse 
(GB)(1989). 


The quantity Nu—1 is a measure of the convective heat 
flux across the fluid layer, where the Nusselt number is 
defined as: Nu = 1 + Pr"^ < V 36 >xyz- The symbol 
<>xyz stands for the spatial average. Figure [5] shows 
the scaling behaviour of the time averaged value of the 
convective heat flux with e for Pr = 0.025. The points 
show the results obtained from the model for different 
values of Q and the dashed line is parallel to the best 
fit over data points near the onset. The best fit shows 
that the convective heat flux scales with e as near 
the primary instability. There is a sharp fall in the heat 
flux at the onset of time-dependent convection. The time 
averaged value of the convective heat flux then starts in¬ 
creasing initially much faster and then much slowly with 
increase in e. The variation of the convective heat flux 
with e becomes almost identical for all values of Q inves¬ 
tigated at higher values of e. In this regime, the slope 
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the Nu— e curve at higher values of e is lower than 0.91. 
The continuous curves show the variation of convective 
heat flux with e, as reported by Clever and Busse [l^ 
for Pr = 0.025 with no-slip boundary conditions. They 
found the scaling exponent to be much larger (> 1.4) for 
the stationary magnetoconvection, which is unusual. The 
qualitative behaviour of the results obtained from the 
model has broad similarity with those obtained by Clever 
and Busse [T^, but our model always shows the scaling 
exponent of heat flux with e less than unity near the onset 
of stationary convection. It is also in excellent agreement 
with the recent results from direct numerical simulations 
with free-slip boundary conditions No experiment 

suggests that the scaling behaviour Nu — 1 ~ e“ with 
the scaling exponent a > 1.4 near the primary instabil¬ 
ity (stationary convection). Like other global variables, 
the heat flux also decreases at the onset of oscillatory in¬ 
stability, reaches a minimum and then increases with e. 
The scaling exponent is less than unity for e > l(r > 2). 



e 


FIG. 7: Plot of convective heat flux {Nu — 1) versus e = 
[Ra — Rac{Q)]/Rac{Q) for Pr = 0.1 as obtained from the 
model (M) and obtained by Clever and Busse (CB)(1989). 




Immediately after the secondary instability, K increases 
linearly with e. Figure [HKb) displays the variation of $ 
with e for Pr = 0.1. The global variable $ scales with 
e as ® near the primary instability. All the curves for 
4)(e) again merge for 1 < e < 10. The scaling exponent 
of 4) in this region is found to be 0.4. 

Figure 0 shows the plot of — 1 versus e for Pr = 0.1. 
The points connected by dots are the results obtained 
from the model for different values of Chandrasekhar’s 
number {Q = 0,30,50 and 100). Three solid lines are 
the results from DNS (Clever and Busse) [l^ with no¬ 
slip conditions for Q = 30, 50 and 100. The linear dashed 
line is the best fit of data points obtained from the model 
for magnetoconvection near the onset. The convective 
heat flux shows the same scaling behaviour as before in 
the regime of stationary magnetoconvection. The results 
obtained from the model are in better agreement with 
those of Clever and Busse [l2| for Pr = 0.1. For higher 
values of e, all the curves seem to become parallel to each 
other. This is in agreement with the results of DNS. 


VI. TIME-DEPENDENT FLUID PATTERNS 


FIG. 6: Variations of (a) the mean kinetic energy K and (b) 
the mean convective entropy 4? with e for Pr = 0.1 as obtained 
from the model. 

Figure [HKa) shows the variation of mean kinetic en¬ 
ergy K as a function e for Pr = 0.1 and for different 
values of Q. The kinetic energy K varies again linearly 
with e near the onset of the stationary convection in the 
form of straight rolls, and decreases at the onset of time- 
periodic instability. However, oscillatory instability oc¬ 
curs at much higher values or r in this case, and the de¬ 
crease in K at the onset of oscillatory instability is small. 


As the Rayleigh number is increased, different time de¬ 
pendent patterns appear at the secondary instability for 
different values of Pr. The threshold for the secondary 
instability increases as Q is raised. For low-Prandtl- 
number convection P [ 2 ^, the secondary instability is 
oscillatory convection, which occurs close to the onset 
of primary convection. This suggests that a qualitative 
description of the convection near secondary instability 
in the presence of small magnetic field may be captured 
qualitatively with relatively less number of modes. The 
selection of the modes for the vertical vorticity are done 
to construct a minimum mode model to be able to cap¬ 
ture the scaling behaviour of the global quantities of 
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Fluid 

patterns 

Pr = 0.025 

Pr = 0.1 

r (Q = 0) 

r (Q = 10) 

r (Q = 30) 

r (Q = 50) 

r (Q = 0) 

r (Q = 10) 

r (Q = 30) 

r (Q = 50) 

2D Rolls 

< 1.04 

< 1.05 

< 1.09 

< 1.18 

< 1.54 

< 1.62 

< 1.82 

< 2.07 

WR 

- 

- 

- 

- 

1.55 - 1.62 

1.63 

- 

- 

QWR 

1.05 - 1.11 

1.06 - 1.12 

1.10 - 1.14 

- 

1.63 - 1.64 

1.64 

- 

- 

CWR 

1.12 - 1.14 

1.13 - 1.27 

1.15 - 1.26 

1.19 - 1.26 

1.65 - 1.75 

1.65 - 1.77 

1.83 - 1.85 

- 

QWR 

1.15 - 1.18 

- 

- 

- 

- 

- 

- 

- 

CWR 

1.19 - 1.29 

- 

- 

- 

- 

- 

- 

- 

TW 

1.30 - 1.31 

1.28 - 1.30 

1.27 - 1.29 

1.27 - 1.29 

1.76 - 1.77 

1.78 - 1.79 

1.86 - 1.88 

- 

QTW 

- 

- 

1.30 

1.30 

1.78 - 1.92 

1.80 - 1.95 

1.89 - 2.04 

2.08 - 2.13 

CTW 

> 1.32 

> 1.31 

> 1.31 

> 1.31 

> 1.93 

> 1.96 

> 2.05 

> 2.14 


TABLE II: Magnetoconvective patterns obtained from the model for Pr = 0.025 and Pr = 0.1. Two-dimensional stationary 
rolls (2D Rolls), periodic wavy rolls (WR), quasiperiodic wavy rolls (QWR), chaotic wavy rolls (CWR), periodic traveling 
waves (TW), quasiperiodic traveling waves (QTW) and chaotic traveling waves (CTW) are observed in the model. 


Rayleigh-Benard magnetoconvection near the onset. The 
scaling properties are captured quite well in the model. 
The fluid patterns are however sensitive to Pr, r, Q and 
wavenumber of the perturbations. A single model is un¬ 
likely to capture the patterns with variations of all these 
parameters. In addition, we have considered only the 
critical wave number kc{Q). We have tested the model 
by adding more vorticity modes. The oscillatory nature 
of the secondary stability does not change qualitatively. 
However, the onset of higher order (e.g. tartiary) insta¬ 
bility is affected. We now discuss the magnetoconvective 
patterns obtained from the model presented here. 

We observe temporal quasiperiodic wavy rolls (QWR) 
for Pr = 0.025 and periodic wavy rolls (WR) for Pr = 
0.1 at the onset of secondary instability. In the latter 
case, the convection becomes temporally quasiperiodic 
with increase in r. Further increase in r leads to chaotic 
wavy rolls (CWR). The onset of chaotic waves also de¬ 
layed for larger value of Q. For a fixed value of Q the 
onset of secondary instability is higher for higher value 
of Pr. At relatively higher values r, the time depen¬ 
dent convective flow bifurcates from standing waves to 
traveling waves. The convective flow consists of periodic 
traveling waves (TW) in a narrow range of r, which de¬ 
pends on Pr and Q. Wavy rolls travel in a direction 
perpendicular to the direction of roll-axis. Further in¬ 
crease in r leads to quasiperiodic traveling waves (QTW) 
and chaotic traveling waves (CTW). Table HIl lists the 
fluid patterns observed for different values of r, Pr and 

Q. 

We now discuss the results on the fluid patterns ob¬ 
served for Pr = 0.1 and Q = 10. Figure [5] shows the 
contour plots of the convective temperature field at mid¬ 
plane z = 0 for r = 1.63 at four different instants. The 
convection shows wavy rolls (WR). They represent stand¬ 
ing wave solutions. The positions of nodes do not vary in 
time. The dimensionless period oscillation for this peri¬ 
odic standing wave solution is r = 0.88. Figure [9] shows 
the mid-plane contour plots of the convective tempera- 
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FIG. 8: Contour plots of the convective temperature field 
at mid-plane 2 = 0 showing wavy rolls (WR) for Pr — 0.1, 
Q = 10 {kc = 3.25) and r = 1.63 at four instants (a) t = 0, 
(b) t = r/4, (c) t = T12.. and (d) i = 3r/4, where r = 0.88 is 
the dimensionless time period of the wavy rolls. 


ture field for r = 1.78 at four instants (a) t = 0, (b) 
t = 0.40, (c) t = 0.80, and (d) t = 1.20. Convective 
patterns are shown to travel slowly along the positive 
X— direction. The inclined arrows mark the position of 
the same hot (white) region which now moves with time. 
The structures along the y-axis make the patterns three 
dimensional, which travel from left to right in this view- 
graph. 

As r is raised further, the periodic traveling waves be¬ 
come quasiperiodic in time. Figure [10] shows the tem¬ 
poral variation of the two largest Fourier modes: (a) 





































FIG. 9: Contour plots of the convective temperature field at 
the mid-plane z = 0 showing periodic traveling waves (TW) 
for Pr = 0.1, Q = 10 (fcc = 3.25) and r = 1.78 at (a) t = 0, (b) 
t = 0.40, (c) t = 0.80, and (d) t = 1.20. The inclined arrows 
show the positions of the same hot (white) region with time. 
The fluid patterns move towards the right. 
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FIG. 11: The projection of the phase space on the Wm — 
Wioi plane showing quasiperiodic magnetoconvection. All 
parameters are same as given in Fig. 1101 


lar region of the phase space gets filled completely. This 
confirms the convection to be quasiperiodic. 

Figure [12] displays the mid-plane contour plots of the 
convective temperature field for this case {r = 1.9) at dif¬ 
ferent instants of time. The patterns consists of inclined 
wavy rolls. The axis of the wavy rolls keep alternat¬ 
ing in time between mutually perpendicular directions. 
The system of alternating oblique wavy rolls also keep 
traveling along the negative direction of the x-axis. The 
inclined arrows shown in the figure follow the positions 
of the same hot region, which moves towards the left as 
time passes. This complex and new magnetoconvective 
patterns are quasiperiodic in time. We have not investi¬ 
gated the possibility of convection in the form of square 
patterns. They may be possible at relatively hi ghe r val¬ 
ues of r. In the DNS with stress-free boundaries [^, two 
sets of stationary rolls leading to patterns of asymmetric 
squares (cross-rolls) appear at relatively higher values r 
at the tartiary instability. 


FIG. 10: The temporal variation of the two largest Fourier 
modes (a) ITioi and (b) Win for Pr = 0.1, Q — 10 and 
r = 1.90 [kc(Q) = 3.25] showing a quasiperiodic. 


Wioi and (b)kFiii for r = 1.9, Pr = 0.1 and Q = 10. 
These signals suggest quasiperiodic magnetoconvection. 
The period of amplitude modulation for the Fourier mode 
Will is almost double of the time of amplitude modula¬ 
tion of the Fourier mode Wim. Figure [TT] shows the pro¬ 
jection of the phase space of the hydro-magnetic system 
on the VFiii — Wioi plane. As time passes, a particu- 


VII. CONVECTIVE TEMPERATURE AND 
VELOCITY PROFILES 

We now present the results on the temperature and 
velocity profiles obtained from the model. We define 
a symbol < ... >x,y,t to describe the spatial average 
in horizontal plane and the temporal average of any 
quantity inside the angular bracket. We then study the 
root mean square (rms) of any relevant quantity / as 
/ (rms) = i/< >x,y,t, which is a function of the ver¬ 

tical coordinate z. Figure [PTI shows the variations of (a) 
vi (rms) and (b) V 3 (rms) along the vertical axis for sta¬ 
tionary straight rolls in fluid with Pr = 0.1 at r = 1.05 
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FIG. 12: Contour plots of the convective temperature field at 
mid-plane 2 = 0 for r = 1.9 showing quasiperiodic traveling 
waves (QTW) at several instants: (a) t = 0, (b) t = 0.60, (c) 
t = 1.20, (d) t = 1.80, (e) t = 2.40, and (f) t = 3.00. The 
inclined arrows show the positions of the same hot region with 
time. Other parameters are as given in Fig. 1101 


for different values of Q. V 2 (rms) is always zero for 
straight rolls parallel to j/-axis. The peak of the hori¬ 
zontal velocity field [vi (rms)] decreases and that for the 
vertical velocity field [ua (rms)] increases with increase 
in Q. The variation in the rms value of the convective 
temperature 6 (rms) in the vertical direction for straight 
rolls is shown in Fig. [Ij] (c). The rms of the convective 
temperature has a peak in the middle of the cell just 
above the onset of convection, as in the case of the verti¬ 
cal velocity. However, the peak rms value decreases with 
increase in Q value. The temporal mean of horizontally 
averaged convective temperature field is more or less si¬ 
nusoidal [Fig. [13] (d)]. The effect of the vertical magnetic 
field is negligible for the Q values investigated here. 

Figure [T3| shows the variations of vi (rms) (blue 
dashed line with points), V 2 (rms) (green dotted line 
with points), and V 3 (rms) (red solid line with points) 
along the vertical direction for Pr = 0.1, Q = 10 
(kc(Q) = 3.25) and for different values of values of r. Fig¬ 
ure [Tdl' a'l describes the case of stationary straight (2D) 
rolls along the y-axis for r = 1.50, where V 2 (rms) is al¬ 




FIG. 13: Variations of (a) v\ (rms) and (b) V3 (rms), (c) 
6 (rms), and (d) < 9 > along 2 — axis for stationary 2D con¬ 
vection (Pr — 0.1, r = 1.05) for different values of Q. The 
velocity V 2 (rms) is always zero for 2D rolls parallel to the 
2 /-axis. 



FIG. 14: Variations of root mean square velocities in the verti¬ 
cal ( 2 ) direction for four different values of r: (a) r = 1.50, (b) 
r = 1.63, (c) r = 1.78, (d) r = 1.90. ui (rms)(z), V 2 (rms)(z) 
and V 3 (rms)(z) are represented by blue dashed-dotted curve, 
green line with dots and red curve with dots, respectively. 
Other parameters are: Pr — 0.1, Q = 10 (fcc = 3.25). 


ways zero. vi (rms) has two peaks: one in the upper 
part and another in the lower part of the cell. It van¬ 
ishes in the middle of the cell (z = 0). The rms value 
of the vertical velocity V 3 (rms) has a peak at z = 0. 
All quantities are symmetric about 2 = 0 and vanish at 
the top and bottom boundaries 2 = ±0.5 due to no-slip 
conditions. Fig. [HI bl describes three-dimensional con- 
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FIG. 15: Variation of root mean square value of the convective 
temperature field in the vertical direction for four different 
values of r: (a) r — 1.50, (b) r = 1.63, (c) r = 1.78, (d) r = 
1.90. Other parameters are: Pr = 0.1, Q = 10, kc{Q) — 3.25. 


vection for r = 1.63. The velocity along the roll-axis 
V 2 {rms) becomes non-zero, as three-dimensional convec¬ 
tion begins. It shows two peaks located at the positions 
of peaks in vi {rms), but with less in magnitude. The 
minima of vi {rms) at 2 = 0 is now non-zero while the 
minima of V 2 {rms) is zero in the middle of the cell. The 
peak of the vertical velocity V 3 {rms) at z = 0 is slightly 
reduced. With further increase in r, the peak values of 
V2 {rms) and V3 {rms) [Figs. nTT cl and (d)] also increase, 
while that of vi {rms) remains almost constant. 

Figure [15] shows the variation of the rms value of the 
convective temperature in vertical direction as r is varied 
keeping all other parameters fixed {Q = 10 and Pr = 0.1) 
for convective structures. The rms of the convective tem¬ 
perature, which had a peak in the middle at the onset of 
stationary convection, becomes flat with increase in r. 
Figure [TS] (a) shows the rms of the convective tempera¬ 
ture at r = 1.50. The rms value value of the convective 
temperature drops in the middle of the cell at the onset 
of oscillatory convection [Fig. [15] (b)], which is expected. 
As r is raised further, the convective temperature shows 


bimodal behavior [Figs.[Tni(c) and (d)]. 


VIII. CONCLUSIONS 


We have presented in this paper a low dimensional 
model for convection of an electrically conducting fluid 
enclosed between two rigid horizontal boundaries in a 
uniform vertical magnetic field. We find that the global 
quantities scale with the distance from criticality e = r — 1 
near the onset of primary instability. The kinetic energy 
K , convective entropy and convective heat flux Nu—1 
scale with e as e, ® and respectively. The vertical 
magnetic field inhibits the primary as well as secondary 
instabilities. Onset is always in the form of 2D station¬ 
ary convection for the range of Pr investigated. The 
secondary instability is time dependent, which may ei¬ 
ther be periodic or quasiperiodic. The model shows the 
appearance of standing waves at the onset of secondary 
instability. With further increase of r, there is a bifur¬ 
cation to traveling wave solutions. A traveling wave so¬ 
lution consists of oblique wavy rolls whose axis alternate 
quasiperiodically in time between two directions perpen¬ 
dicular to each other. We find qualitative similarity of 
the results of this simple model with the results of DNS 
by Clever and Busse [l^- The rms value of the verti¬ 
cal velocity has one peak in the middle of the cell. The 
rms values for the horizontal velocities show bi-modal be¬ 
haviour with two peaks are located symmetrically about 
the mid-plane (z = 0). The rms value of the convective 
temperature is flat near the middle of the cell for station¬ 
ary convection and become bi-modal for time dependent 
convection. The two peaks of the convective temperature 
are also located symmetric about the mid-plane. How¬ 
ever, the locations of the peaks for the convective tem¬ 
perature is closer to the mid-plane compared with the 
same for the peaks of the horizontal velocities. The rms 
values of all fields vary with variation of r and Q. 
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